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Abstract 

We investigate baryon-baryon (BB) interactions in the 3-flavor full QCD simula- 
tions with degenerate quark masses for all flavors. The BB potentials in the orbital 
S-wave are extracted from the Nambu-Bethe-Salpeter wave functions measured on the 
lattice. We observe strong flavor-spin dependences of the BB potentials at short dis- 
tances. In particular, a strong repulsive core exists in the flavor-octet and spin-singlet 
channel (the 8 S representation), while an attractive core appears in the flavor singlet 
channel (the 1 representation). We discuss a relation of such flavor-spin dependence 
with the Pauli exclusion principle in the quark level. Possible existence of an H- 
dibaryon resonance above the AA threshold is also discussed. 



§1. Introduction 

The generalized nuclear force, which includes not only the nucleon-nucleon (NN) interac- 
tion but also hyperon-nucleon (YN) and hyperon-hyperon (YY) interactions, has been one 
of the central topics in hadron and nuclear physics. Experimental studies on the ordinary 
and hyper nuclei 1 - 1 as well as the observational studies of the neutron stars and supernova 
explosions 2 -* are intimately related to the physics of the generalized nuclear force. 

As far as the NN interaction is concerned, significant number of scattering data have been 
accumulated. Furthermore, the phenomenological NN potential is one of the useful means to 
parametrize the phase shift data below the meson production threshold: Indeed, the S-wave 
NN phase shift can be described well by a potential in the coordinate space with a repulsive 
core at short distance and an attractive well at intermediate and long distances. 3 - ) ~ 5 - ) 

In contrast to the NN interaction, properties of the YN and YY interactions are poorly 
known even today due to the lack of high precision scattering data for hyperons. The 
approximate flavor SU(3) symmetry in the hadronic level does not fully constrain the hyperon 
interactions since there exists six independent flavor-channels for the scattering of octet 
baryons. Theoretical attempts to extract the YN and YY interactions have been made on 
the basis of the constituent quark models 6 -* or the one-boson-exchange models, 7 - 1 but there 
still exists a large uncertainty for the determination of the generalized nuclear force. To 
improve this situation, a model independent investigation based on QCD is called for. 8 ) 

Recently a new method to extract the baryon-baryon potentials in the coordinate space 
from lattice QCD simulations has been proposed and applied to the NN system, 9 )' 10 ) and EN 
and AN systems. 11 - 1 ' 12 > In this method, energy independent non-local potentials are defined 
from the Nambu-Bethe-Salpeter (NBS) amplitude through the Schrodinger equation. In the 
above studies, it has been shown that the repulsive core at short distance appears in all NN, 
EN and AN systems, while the height and range of the repulsion depend on flavor and spin 
combinations of two baryons. 

The purpose of this paper is to make a systematic study on the short range baryon-baryon 
(BB) interactions with full QCD simulations on the lattice. (A preliminary account is given 
in Ref. 13 )). This problem is also closely related to possible dibaryon states such as the flavor- 
singlet six-quark system H in phenomenological quark models. 14 - 1 ' 15 -* To capture essential 
features of the short range force without contamination from the quark mass difference, we 
consider the flavor SU(3) limit where all u, d, and s quarks have a common finite mass. This 
allows us to extract BB potentials for irreducible flavor multiplets, from recombinations of 
which the potentials between asymptotic baryon states are obtained with suitable Clebsch- 
Gordan (CG) coefficients. 



This paper is organized as follows. In Section 2, we explain a classification of BB states 
in the flavor SU(3) limit and a method to extract BB potentials in lattice QCD. In Section 3, 
we give a brief summary of our numerical simulations. In Section 4, we present the resulting 
BB potentials and discuss their implications. The last section is devoted to summary and 
outlook. 

§2. Baryon-baryon potentials in the flavor SU(3) limit 

In the flavor SU(3) limit, the ground state baryons belong to the flavor-octet with spin 
1/2 or the flavor-decuplet with spin 3/2. In this paper, we focus on the octet baryons, for 
which two baryon states with a given angular momentum are labeled by the irreducible flavor 
multiplets as 

8 <g> 8 = 27 © 8 S © 1 © 10* © 10 © 8 a . (1) 

symmetric anti-symmetric 

Here "symmetric" and "anti-symmetric" stand for the symmetry under the flavor exchange 
of two baryons. For the system in the orbital S-wave, the Pauli principle between two baryons 
imposes 27, 8 S and 1 to be spin singlet ( 1 5'o) while 10*, 10 and 8 a to be spin triplet ( 3 Si). 
Since there are no mixings among different multiplets in the SU(3) limit, one can define the 
corresponding potentials as 

% : \/ (27) (r), V {8 °\r), V w (r), (2) 

3 ^ : V (1(r) (r), V {w \r), V {8a \r) . (3) 

Potentials among octet baryons, both the diagonal part (BiB 2 — > BiB 2 ) and the off- 
diagonal part (BiB 2 — > B 3 fi 4 ), are obtained by suitable combinations of V^ a \r) with 
a = 27,8 s ,l,10*,10,8 a . 

According to Ref. 9 -'' 10 -' , the above potentials can be defined from the NBS wave function 
<p( a \r) through the Schrodinger equation as 



H(f)= /dW (a) (f,O0 (a) (O> (4) 



\7 2 

— + E™ 

2/i 

with E^ = k 2 /(2fi) and k 2 = W 2 /A - M 2 /A. Here // = m/2, M = 2m and W are the 
reduced mass, total mass and the total energy in the center of mass frame of the two baryon 
system with m being the octet baryon mass in the SU(3) limit. Note that U^(f,f) is non- 
local but energy-independent and may be written as U^(r, f") = V^(r, V)5 3 (r — -P). 10 ) 
At the leading-order of the derivative expansion of V^(f, V), 10 -*' 16 - 1 we obtain an effective 



Table I. Summary of lattice parameters and hadron masses. The value of lattice spacing a is 
determined from p meson mass at the physical light and strange quark masses. See the official 
home-page of CP-PACS and JLQCD Collaborations. 18 ) 
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(0\(BB)W(t,r)\B = 2,a,W), 



(5) 



(6) 



where (BB)( a \t,r) = ^2 i • -CijBi(t,x + r)Bj(t,x) is a two-baryon operator with relative 
distance r*in a-plet and \B = 2, a, W) is a QCD eigenstate with the total energy W and the 
baryon number 2 in the ct-representation. Here B t is an one-baryon composite field operator 
in the octet. The relation between two-baryon operators in the flavor basis and baryon basis 
are given in Appendix A. In the lattice QCD simulations, the above NBS wave function for 
the smallest W is extracted from the four point function as 



G { ±\t -t ,r) = (0| (BBY a \t,r) (BB) [a \t ) |0) oc ^ a \r) e - w ^- to) 



(7) 



r("), 



for t — to ^> 1, where (BB) (t ) is a wall source operator at time to to create two-baryon 
states in a-plet while (BB)^(t,f) is the sink operator at time t to annihilate two-baryon 
states. 



§3. Numerical simulations 

We employ the gauge configurations generated by CP-PACS and JLQCD Collaborations 
with the renormalization group improved Iwasaki gauge action and the non-perturbatively 
O(a) improved Wilson quark action. 17 ) These configurations are provided by Japan Lattice 
Data Grid (JLDG) and International Lattice Data Grid (ILDG). Quark propagators are 
calculated for the spatial wall source at to with the Dirichlet boundary condition in temporal 
direction at t — to ± 16 mod 32. In order to enhance the signal, all 32 time slices on each 



*' The "effective" central potential is obtained from the NBS wave function without making a decom- 
position into central and tensor potentials in the leading order of the derivative expansion. 9 '' 10 ' In the 
^o state, the "effective" central potential reduces to the central potential due to the absence of the tensor 
component. 
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Fig. 1. NBS wave functions at m^ = 835 MeV and t — t$ = 10, normalized to 1/2 for the singlet 
channel and to 1 for other channels at the maximum distance. 



configuration are used as to of the wall source. In addition we carry out an average over 
forward and backward propagations in time. 

Our lattice parameters are summarized in Table I. Those hopping parameters reside 
around the physical strange quark mass region in 2+1 flavor QCD at the same /3. 18 ^ Since 
the typical range of the BB interactions at this quark mass is about 1 fm or less as shown 
below, we expect that the small lattice size (L ~ 2 fm) would not affect the short range part 
of the BB interactions qualitatively. 

In our calculation, the sink operator is projected to the Af representation of the cubic 
group, so that the NBS wave function is assumed to be dominated by the S-wave component. 
We estimate the statistical errors by the jackknife method with a bin size of seven for 
K u ds = 0.13710 and eight for K u d s = 0.13760. All of these numerical computations have been 
carried out at KEK supercomputer system, Blue Gene/L and SR11000. 

§4. Results and their implications 

4.1. BB potentials in flavor basis 

Fig. 1 shows the NBS wave functions as a function of the relative distance between 
two baryons at m^ = 835 MeV and t — to = 10. To draw all data in a same scale, they 
are normalized to 1/2 for the singlet channel and to 1 for other channels at the maximum 
distance. Since the effective mass of the four point function in each channel shows a plateau 
at t — to > 10, we use data at t — to — 10 exclusively to extract both NBS wave functions 
and BB potentials throughout this paper. The wave functions in Fig. 1 show characteristic 
flavor dependence: In particular, a strong suppression at short distance appears in the 8 S 
channel, while a strong enhancement appears in the 1 channel. Similar results are obtained 
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2. The six independent BB potentials for S-wave in the flavor SU(3) limit, extracted from the 
lattice QCD simulation at m„- = 1014 MeV (red bars) and m n = 835 MeV (green crosses). 



for m w = 1014 MeV too. 

Fig. 2 shows the resulting six BB potentials in the flavor basis obtained from the NBS 
wave functions. Red bars (green crosses) data correspond to the pion mass 1014 MeV (835 
MeV): Although there is a tendency that the magnitude (range) of the potentials becomes 
larger at short distances (longer at large distances) for lighter quark mass, the differences 
are not substantial for the present heavy quark masses. 



Table II. Baryon pairs in an irreducible flavor SU(3) representation, where {BB'} and [BB'\ 
denotes BB' + B'B and BB' - B'B, respectively. 



flavor multiplet 


baryon pair (isospin) 


27 


{NN}(1=1), {NZ}(I=3/2), {EE}{I=2), 




{£~}(I=3/2), {SS}(I=1) 


8 S 


none 


1 


none 


10* 


[NN] (1=0), [E~] (1=3/2) 


10 


[NE] (1=3/2), [~ ~](I=0) 


8 a 


[TVS] (1=0) 



The left three panels show potentials in the 1 So channel, while the right three panels 
show effective central potentials in the 3 Si channel. We have estimated the energies E^ 
from the condition that Eq.(5) approaches zero at large r (the method called "from V" in 
Ref. 19 )). We find E^ ~ -5 MeV, E^ ~ 25 MeV, #W ~ -30 MeV, £( 10 *) ~ -10 MeV, 
_g( 10 ) r^ q M e v and E^ 8 ^ ~ — 15 MeV. Since the constant term (the second term in the right 
hand side of Eq.(5)) is small compared to the Laplacian term (the first term in the right 
hand side of Eq.(5)) at short distance, we have not attempted to extract a precise value of 
E^ in the present study. 

Some of octet-baryon pairs belong exclusively to one irreducible representation, as shown 
in Table II. For example, a symmetric NN belongs to 27 so that one can regard \/( 27 ) as 
a SU(3) limit of the NN 1 5'o potential. Similarly, V^ 10 *\ V^ 10 ^ and V^ 8 "-* can be considered 
as the SU(3) limits of some potentials in the baryon basis, while V^ 8 ^ and V^ are always 
superpositions of different potentials in the baryon basis. 

Top two panels of Fig. 2 show \A 2T ) and \A 10 *), which corresponds to NN 1 5'o potential 
and NN 3 S± potential respectively. Both have a repulsive core at short distance and an 
attractive pocket around 0.6 fm. These qualitative features are consistent with the previous 
results found in the NN system in quenched approximation with lighter quark mass. 9 - 1 The 
right middle panel of Fig. 2 shows that V^ 10 ^ has a stronger repulsive core and a weaker 
attractive pocket than V^ 27 ' 10 *- 1 . Furthermore V^ 8 ^ in the left middle panel of Fig.2 has 
a very strong repulsive core among all channels, while V^ 8 ^ in the right bottom panel of 
Fig.2 has a very weak repulsive core. In contrast to other cases, V^ shows attraction for all 
distances instead of repulsion at short distance, as shown in the left bottom panel of Fig. 2. 

Above features are consistent with what has been observed in phenomenological quark 
models. 6 -* In particular, the potential in the 8 S channel in quark models becomes strongly 



repulsive at short distance since the six quarks cannot occupy the same orbital state due to 
quark Pauli blocking. On the other hand, the potential in the 1 channel does not suffer from 
the quark Pauli blocking and can become attractive due to short range gluon exchange. Such 
an agreement between the lattice data and the phenomenological models suggests that the 
quark Pauli blocking plays an essential role for the repulsive core in BB systems as originally 
proposed in Ref. 20 -* 

4.2. BB potentials in baryon basis 

The BB potentials in the baryon basis can be obtained by an unitary rotation of those 
in the flavor basis. Its explicit form reads, 

a 

where U is an unitary matrix which rotates the flavor basis {{a\} to baryon basis {{i\}, i.e. 
(i\ = Y^ a Ui a {a\. The explicit forms of the unitary matrix U in terms of the CG coefficients 
are given in Appendix B. 

In Fig. 3, as characteristic examples, we show the potentials for S=— 2, 1=0 sector in the 
1 S'o channel at m n = 835 MeV. To obtain V^-(r), we first fit the potentials in the flavor basis 
by the following form with five parameters b\ ~ b 5 , 

/ ^ ;, 5 r \ 2 

V(r) = b 1 e- b2r2 + b 3 (l-e- b ^ 2 )(^—) . (9) 



We then use the right hand side of Eq.(8) to obtain the potentials in the baryon basis. The 
left panel of Fig. 3 shows the diagonal part of the potentials. The strong repulsion in the 8 S 
channel is reflected most in the EE(l=0) potential due to its largest CG coefficient among 
three channels. The strong attraction in the 1 channel is reflected most in the NS(l=0) 
potential due to its largest CG coefficient. Nevertheless, all three diagonal potentials have 
repulsive core originating from the 8 S component. The right panel of Fig. 3 shows the 
off-diagonal part of the potentials which are comparable in magnitude to the diagonal ones. 
Since the off-diagonal parts are not negligible in the baryon basis, full coupled channel 
analysis is necessary to study observables. Similar situation holds in (2+l)-flavors: The 
flavor basis with approximately diagonal potentials are useful for obtaining essential features 
of the BB interactions, while the baryon basis with substantial magnitude of the off-diagonal 
potentials are necessary for practical applications. 

Other potentials in baryon basis are given in Fig. 4. Since the 8 S state does not couple 
to the 3 Si channel, the repulsive cores in the 3 S\ channel (the right panels of Fig. 4) are 
relatively small. The off-diagonal potentials are not generally small: In the right second panel 
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Fig. 3. BB potentials in baryon basis for S=— 2, 1=0, 1 Sq sector. Three diagonal(off-diagonal) 
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of Fig. 4, for example, the NA-NE potential in the 3 Si channel is comparable in magnitude 
at short distances with the diagonal NA-NA and NE-NE potentials. Although all quark 
masses of 3 flavors are degenerate and rather heavy in our simulations, these coupled channel 
potentials in the baryon basis may give useful hints for the behavior of hyperons (A, E and 
S) in hypernuclei and in neutron stars. 1 '' 2 ' 

§5. Summary and outlook 

We have performed 3-flavor full QCD simulations to study the general features of the 
BB interaction in the flavor SU(3) limit. From the NBS wave function measured on the 
lattice, we extracted all six independent potentials in the S-wave at the leading order of 
the derivative expansion; V^ 27 \ V^ 8s \ V^\ V^ 10 *\ V^ 10 ^ and V^ labeled by flavor irre- 
ducible representation. We have found strong flavor dependence of the BB interactions. In 
particular, V^ 8 ^ has a very strong repulsive core at short distance, while V^ 1 ' is attractive 
in all distances. These features are qualitatively consistent with the Pauli blocking effect 
among quarks previously studied in phenomenological quark models. Recently, the existence 
(absence) of the repulsive core in the K + -P potential 21 ) (in the cc-N potential 22 ^) has been 
observed in the same method as adopted in this paper. These results further support the 
relevance of the Pauli blocking effect in the quark level for hadron interactions. 

The BB potentials in the baryon basis are reconstructed from those in the flavor basis. 
We have found that the potentials in 3 5*i channels have weaker repulsive core than the 
1 So channels, since the 8 S state, which has the strongest repulsive core, does not mix with 
the former. We note that the off-diagonal potentials are not generally small compared 
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Fig. 4. BB potentials in baryon basis other than those shown in Fig. 3. See the caption of Fig. 3. 
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to the diagonal ones. These information together with the future lattice simulations at 
non-degenerate quark masses in larger volumes will give useful hints for the structure of 
hypernuclei and neutron stars. 

The flavor singlet channel has attraction for all distances in our simulation. The present 
data are not sufficient to make a definite conclusion on the if-dibaryon, because of single 
small lattice volume and a large degenerate quark mass. Nevertheless, by using the potential 
V^(r) and solving the Schrodinger equation Eq.(4) with E^ = —30 MeV, we find a shallow 
bound state, which suggests a possibility of a bound iJ-dibaryon with baryon components 
AA : EE : SN = — vT : v^ : \fi in the SU(3) symmetric world at large quark masses. In 
the real world, the baryon mass ordering in S=— 2 and 1=0 sector becomes AA < SN < EE 
due to flavor-SU(3) breaking. In our trial analysis of solving the Schrodinger equation by 
using the potentials obtained in this paper together with small baryon mass differences from 
a 2+1 flavor lattice QCD simulation*), a resonance state is found at an energy between AA 
mass and EN mass. If such a resonance exists in nature, it may explain the enhancement 
just above the AA threshold recently reported in the KEK experiment. 23 ) However, further 
investigations in both theory and experiment are necessarily to draw a definite conclusion. 
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Appendix A 

The irreducible BB operators 

In this Appendix, flavor irreducible BB operators used in this study are given. The 
composite operators for octet baryons are 



Pa(x) = +e cl , C2 , ca (C , 75) j 8i,a,<Wu(£lM6)u(&) 


(A-l) 


n a ( X ) = +e Cl , C2 , C3 (<^75)/3 1 ,/3 2 ^ 3 , a M(6)^(6)^fe) 


(A-2) 


j 


(A-3) 


-£0*0 = - e ci,c 2 ,c 3 (C75W 2 <W y- K6)s(6)m(6) + M6M6M&)] 


(A-4) 


-^fa) = -W2.es (^Wft W rf (£lM6M£ 3 ) 


(A-5) 


S a(^) = + e C1 , C2 , C3 (CTs)^,^ fys,a s(6)m(6)s(Cs) 


(A-6) 


^W = +e Cl , C2 ,c3 (C75)/3i,/3 2 <W s(6)^(6)s(6) 


(A-7) 


A a (x) = — €ci,c2,c 3 \Clb)pi,fa $133,01 X 





i K6M6M6) + *(6M6M6) - 2 M (6)rf(6)s(6)] (A-8) 

with a notation £j = {cj,/3j,x}. We follow the phase convention in Ref. 25 ^ . In this con- 
vention, the two-baryon operator which belongs to the definite flavor representation can be 
constructed with the Clebsch-Gordan coefficient of SU(3) as 

/07 / 1 /19 

BB^ = +J—AA-J—ES + J—NS (A-9) 

V 40 V 40 V 40 v ' 

BBM = -<JIaA-JIlE+JInB (A.10) 

bbw = -J±aa+^ze + J1ne (A-ii) 



BB (wt) = +J-pn- J-np (A-12) 

bbW = +J\pE + -J\e+ p (A-i3) 



2" V 2 

i p " vr ' vi" v r 



££( 8 «) = +J-pS~ - J-~-p-J-n~ +J-~°n (A- 14) 



where 17 Z 1 and iV^ represents 



^ = +^/^ + IT - y ^° ^° + \j\z- Z + (A-15) 

NB = +^B- + ^S-p-J]n^-Jl^n. (A-16) 
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Appendix B 

Relation between the flavor basis and baryon basis 



Unitary matrices U which rotates the flavor basis {(a|} to baryon basis {{i\}, are given 
as follows. 



1. S=-l, 1=1/2, l S sector. 



(NA\ 
(NE\ 



2. S=-l, 1=1/2, 3 #i sector. 



(NA\ 
(NE\ 



(27| 

(8,1 



(10* 

(8a I 



3. S=— 2, 1=0, 1 Sn sector. 



', ^o 



/ (AA\ \ 

(2727| 
V (NS\ J 



( 



40 

1_ 

40 

12 

10 



40 



'20 

40 / 



/ (27| \ 

(8 S | 



4. S=-2, 1=1, 1 S f) sector. 



j, i j., 



(JVS| 
(27/l| 



(27| 
(8,1 



5. S=-2, 1=1, 3 Si sector. 



/ {NS\ \ 


/ 


(27/l| 


= 


V (EZ\ ) 


\ 



1 \ 



6 / 



/ (10*| \ 
(10| 

V (8a| / 



6. S=-3, 1=1/2, % sector. 



(AS\ 

(27~| 



7. S=-3, 1=1/2, 3 Si sector. 



(AS\ 

(27-| 



(27| 
(8,1 



(10| 
(8 J 



(B-l) 



(B-2) 



(B-3) 



(B-4) 



(B-5) 



(B-6) 



(B-7) 
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